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A  theoretical  study  of  the  structural,  elastic  and  electronic  properties  of  spinel  LiTi204  anode  has  been 
performed  by  density  functional  theory  (DFT)  plane-wave  pseudopotential  method.  The  independent 
elastic  constants,  shear  modulus  (G),  bulk  modulus  (B),  and  Young’s  modulus  (£)  are  evaluated,  respec¬ 
tively.  The  results  suggest  that  cubic  LiTi204  is  mechanically  stable.  The  G/B  ratio  of  0.584  indicates  the 
ductility  of  LiTi204  is  good.  The  electron  density  difference  of  LiTi204  shows  that  the  02p  orbits  overlap 
effectively  with  Ti3d  ones,  confirming  the  formations  of  strong  covalent  bonds  between  them,  while  Li  is 
fully  ionized  in  the  lattice.  The  formation  enthalpy  for  LiTi204  is  calculated  to  be  -2070.723  ±1.6  kj  mol-1 . 
The  strong  covalent  bonds  between  O  and  Ti  atoms  are  not  only  responsible  for  the  excellent  mechanical 
stabilities  but  also  very  crucial  for  the  thermodynamic  stability  of  LiTi204  compound.  Furthermore,  in 
Li2Ti204  compound,  the  full  occupation  of  16(c)  sites  by  Li+  not  only  leads  to  a  smaller  Ci2  value  but 
also  leads  to  a  much  larger  C44  one.  Therefore,  the  plasticity  and  ductility  of  the  Li2Ti204  become  poor  in 
comparison  to  LiTi204,  while  the  thermodynamic  stability  of  Li2Ti204  can  be  further  improved  after  the 
Li+  intercalation  of  LiTi204. 

©  2011  Elsevier  B.V.  All  rights  reserved. 


1.  Introduction 

There  is  a  strong  and  ever-growing  demand  for  renewable 
energy  resources,  reliable  energy  storage,  and  conversion  devices  in 
order  to  replace  the  traditional  fossil  fuels.  Rechargeable  batteries 
are  considered  as  the  most  successful  tried  and  true  technologies 
that  can  repeatedly  generate  clean  electricity  from  stored  mate¬ 
rials  and  convert  reversely  electric  energy  into  chemical  energy 
[1  ].  Undoubtedly,  lithium-ion  batteries  (LIBs)  are  the  most  popular 
rechargeable  batteries  because  of  the  efficiency,  lightweight,  and 
high  energy  density.  One  of  the  key  safety  issues  in  LIBs  as  power 
battery  would  be  the  dendritic  lithium  growth  on  the  anode  sur¬ 
face  at  high  charging  current  because  the  conventional  carbonous 
materials  approach  almost  0  V  versus  Li/Li+  at  the  end  of  Li  insertion 
[2].  Developing  anode  materials  with  high  safety  is  one  of  the  key 
challenges  for  lithium-ion  batteries.  Spinel  lithium  titanium  oxides 
have  been  considered  as  one  of  the  promising  alternatives  for  the 
commercial  carbon-based  anode  materials  due  to  its  unique  prop¬ 
erties.  The  poor  electric  conductivity  of  Li4Ti50i2  limits  its  high 
rate  application,  but  LiTi204  anode  materials  have  advantages  in 
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electric  conductivity  facilitating  high  rate  performance  [3,4]. 
Recently,  rapid  advances  in  computing  technology  have  made  it 
possible  to  perform  predictions  of  the  properties  of  materials  based 
on  density  functional  theory  (DFT)  [5].  The  electronic  structure 
and  crystal  parameters  of  LiTi204  has  been  reported  according  to 
the  first  principle  [6].  W.  Ra  et  al.  [7]  reported  that  the  electronic 
structural  changes  of  electrochemically  reacted  LiTi204  by  first 
principles  calculations.  L.  Benco  et  al.  [8]  reported  the  average  inter¬ 
calation  voltages  of  LiTi204  as  cathode  materials  for  lithium-ion 
batteries  calculated  from  first  principles  using  the  LAPW  method. 
Besides  the  electronic  properties,  the  elastic  properties  are  also  of 
particular  interest  as  they  determine  the  mechanical  response  of 
a  material  to  external  stress  and  the  macroscopic  mechanical  sta¬ 
bility.  To  our  knowledge,  no  theoretical  study  on  the  mechanical 
stability  of  LiTi204  has  been  reported.  In  this  paper,  it  is  our  tar¬ 
get  to  contribute  understanding  towards  the  relationship  among 
electronic  structure,  thermodynamic  and  mechanical  stabilities. 

2.  Theory  and  computation  details 

The  present  calculations  were  performed  in  the  density  func¬ 
tional  theory  (DFT)  framework  [9]  implemented  in  the  CASTEP 
package  [10].  The  exchange  and  correlation  effects  are  treated 
with  the  Perdew-Wang  functional  in  the  generalized  gradi¬ 
ent  approximation  scheme  (GGA-PW91)  [11].  All  computations 
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Fig.  1.  Perspective  view  for  LiT^C^  crystal:  (a)  primitive  cell  (2  LiT^C^  units)  and 
(b)  conventional  cell  (8  LiTi204  units).  The  small  green  spheres  represent  the  16(c) 
sites.  Lii+xTi204  can  be  formed  when  foreign  and  8(a)  site  lithium  ion  transfer  to 
16(c)  sites.  For  L12T12O4,  all  16(c)  sites  were  occupied  by  Li+.  (For  interpretation  of 
the  references  to  color  in  this  figure  legend,  the  reader  is  referred  to  the  web  version 
of  the  article.) 


were  performed  with  a  plane-wave  energy  cutoff  of  450  eV,  and 
the  sampling  over  Brillouin  zone  was  treated  by  a  (5x5x5) 
Monkhorst-Pack  mesh.  This  set  of  parameters  ensures  that  the 
total  energy  of  system  can  be  accurately  evaluated,  and  the  energy 
convergence  is  within  5.0  x  1 0-7  eV  atom-1 .  The  primitive  and  con¬ 
ventional  cells  for  LiTi204  crystal  are  shown  in  Fig.  1,  in  which 
the  16(c)  sites  and  the  Li+  transfer  direction  were  also  depicted. 
In  the  following  calculations,  the  periodic  boundary  condition  is 
employed. 


3.  Results  and  discussion 


Mechanical  stability  of  materials  mainly  reflects  the  struc¬ 
tural  stability  under  a  particular  deformation  condition.  According 
to  DFT  calculations  combined  with  a  linear  elastic  theory,  all 
mechanical  parameters  of  single  crystal  materials  can  be  calculated 
accurately.  Elastic  constants  of  a  solid,  which  reflect  the  relation¬ 
ship  between  stress  and  strain,  are  very  important  because  they  are 
related  to  various  fundamental  solid-state  properties  such  as  inter¬ 
atomic  potentials,  equation  of  states  (EoS),  and  phonon  spectra  [12]. 
It  is  well  known  that  a  cubic  crystal  has  only  three  independent  elas¬ 
tic  constants  Cu,  C12  and  C44.  Therefore,  a  set  of  three  equations  is 
required  to  determine  all  the  elastic  constants  [13]. 

The  calculated  elastic  constants  (Cn,  Ci2  and  C44)  are  listed  in 
Table  1.  According  to  the  elastic  constants,  the  shear  modulus  (G) 
[14],  bulk  modulus  (B)  [14],  and  Young’s  modulus  (E)  [15,16]  can 
be  derived  by  the  following  expressions: 

B  =  Cu  +32Cu  (1) 


G  = 


Cu 


Ci  2  +  3C44 


(2) 


n  (Cn-C12)(C11+2C12)  (  ^ 

E- - C^l -  (3> 

For  a  cubic  crystal,  the  mechanical  stability  requires  that  all  the 
independent  elastic  constants  should  satisfy  the  following  criteri- 


ons  [17-20]: 

C11+2C12>0  (4) 

C44  >  0  (5) 

Cn-C12>0  (6) 

Cu  >  0  (7) 


Furthermore,  the  bulk  modulus  B  should  also  satisfy  an  additional 
criterion  [17-20]: 

Q 2  <  B  <  Cu  (8) 

The  results  clearly  suggested  that  under  elastic  strain  perturba¬ 
tions  cubic  LiTi204  is  mechanically  stable.  To  further  elucidate  the 
mechanical  properties,  the  G/B  ratio,  which  is  now  widely  accepted 
as  a  predictor  in  first-principles  materials  design,  is  introduced. 
Pugh  [21  ]  concluded  that  G/B  ratio  reflects  the  ductility  of  the  mate¬ 
rials.  The  smaller  the  G/B  ratio  is,  the  more  ductile  the  materials  are. 
From  the  calculated  G  and  B  values,  the  G/B  ratio  is  calculated  to 
be  0.584.  The  critical  value  which  distinguishes  ductile  and  brittle 
characteristics  was  about  0.571  [22].  In  our  work,  the  G/B  ratio  of 
0.584  is  rather  close  to  the  critical  value,  indicating  the  ductility  of 
cubic  LiTi204  is  good.  Besides  the  G/B  ratio,  it  was  recently  reported 
that  the  (Cn-Ci2)  also  plays  an  important  role  in  the  mechanical 
properties  of  materials  [23].  The  smaller  the  (Cn-Ci2)  is,  the  lower 
the  Young’s  modulus  is,  and  the  better  the  plasticity  is.  In  present 
work,  the  calculated  (Cn-Ci2)  is  117.21  GPa,  which  indicate  the 
Young’s  modulus  is  relatively  high  and  the  plasticity  of  the  LiTi204 
is  relatively  poor.  As  the  mechanical  properties  are  representa¬ 
tions  of  microscopic  bondings  of  systems,  to  bridge  the  relationship 
between  them,  the  electronic  structures  were  discussed  in  detail. 

Fig.  2  depicts  the  electron  density  difference  diagram  for  cubic 
LiTi204.  Due  to  the  formation  of  covalent  bonds,  the  electron  dis¬ 
tribution  is  changed,  and  the  positive  (in  blue)  or  negative  (in 
red)  regions  indicate,  respectively  where  the  electron  density  is 
enriched  or  depleted.  From  Fig.  2,  it  can  be  identified  that  cova¬ 
lent  bonds  between  Ti  and  O  atoms  are  formed,  leading  to  some 
electron  transformation  from  Ti  to  O  atoms.  And  lithium  exists  in 
the  crystal  lattice  as  pure  ions.  According  to  the  electron  redistri¬ 
bution  characteristics,  it  can  also  be  confirmed  that  O  atoms  take 
sp3  hybridizations,  and  the  hybrid  orbits  overlap  effectively  with 
the  3d  ones  of  three  nearest  neighbor  Ti  atoms.  The  formations 
of  Ti-0  covalent  bonds  are  very  crucial  for  the  aforementioned 
mechanical  properties.  As  the  Ti-0  bonds  are  evenly  distributed 
in  three  dimensional  framework,  the  mechanical  properties  thus 
show  isotropic  natures  (i.e.  Ex  =  Ey  =  Ez).  Moreover,  the  strength  of 
Ti-0  bonds  is  expected  to  determine  the  resistance  of  materials 
against  uniaxial  tensions  in  the  three  directions.  The  stronger  the 
bond  is,  the  larger  the  Young’s  modulus  is.  The  data  listed  in  Table  1 
clearly  suggested  that  Ti-0  bonds  are  rather  strong  and  the  Young’s 
modulus  possesses  a  large  value  (152.08  GPa).  Besides  the  Young’s 
modulus,  the  bulk  modulus,  which  reflects  the  resistance  of  a  mate¬ 
rial  against  volume  change  under  hydrostatic  pressure,  is  also  very 
large  (125.11  GPa).  On  the  basis  of  Eq.  (1),  it  can  be  figured  out 
that  Cu  plays  an  important  role  in  determining  the  bulk  modulus. 
Therefore,  the  bonding  between  Ti  and  O  atoms  is  also  impor¬ 
tant  for  the  bulk  modulus.  Although  the  Ti-0  bonds  are  rather 
strong,  the  C44  (C44  =  C55  =  C66)  value  is  however  a  little  bit  small. 
In  comparison  to  the  Young’s  modulus,  the  shear  deformation  is 
smaller.  Hence,  the  theoretical  shear  stress  will  be  reached  before 
the  stress  of  crack  tip  exceeds  the  theoretical  tensile  stress.  This 
can  be  deduced  that  the  (1  0  0 )  (0  0 1 ),  ( 0 1  0 )  (0  0 1 ),  ( 1  0  0 )  (0 1  0), 
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Table  1 

Calculated  elastic  constants  and  modulus  (all  in  GPa)  of  Lii+xTi204  (x  =  0, 1). 


Cn 

Cl2 

C44 

11 

& 

Ez 

B 

G 

G/B 

LiTi204a 

203.25 

86.04 

82.78 

152.08 

152.08 

125.11 

73.11 

0.584 

Li2Ti204b 

204.22 

69.29 

104.30 

169.11 

169.11 

114.27 

89.57 

0.784 

a  Li+  occupies  the  8(a)  sites.  The  optimized  lattice  constant  is  8.489  A.  The  error  for  modulus  is  restricted  to  ±4.854  GPa. 
b  Li+  occupies  all  the  16(c)  sites.  The  lattice  constant  after  geometry  optimization  is  8.412  A.  The  error  for  modulus  is  restricted  to  ±4.473  GPa. 


(00 1 )  (0 1  0),  (0 1  0)  (1  00),  (00 1 )  (1  0  0)  slip  systems  maybe  very 
active  in  LiTi204  compounds.  The  slips  of  ( 1  0  0),  (0 1  0),  and  (001) 
crystal  plane  are  expected  to  be  prior  to  the  cleavage  process  of 
the  Ti-0  bonds,  and  it  may  be  the  origin  for  the  good  ductility  of 
LiTi204. 

Moreover,  it  should  be  emphasized  that  there  is  an  intercalation 
process  of  Li+  into  the  spinel  host  structures  during  realistic  applica¬ 
tions,  and  the  mechanical  properties  should  reflect  such  structural 
changes.  According  to  the  previous  report  [24],  the  foreign  and 
8(a)  site  Li+  will  transfer  to  the  16(c)  sites  (Fig.  1(b)),  and  finally 
all  16(c)  sites  were  occupied  by  Li+  to  achieve  a  Li2Ti204  formula. 
Table  1  also  lists  the  corresponding  elastic  constants  and  moduli 
of  Li2Ti204.  Based  on  the  criterion  aforementioned,  Li2Ti204  com¬ 
pound  is  mechanical  stable.  In  comparison  to  LiTi204,  the  C44  value 
increases  to  a  very  large  extent,  and  this  means  that  the  occupa¬ 
tion  of  1 6(c)  sites  by  Li+  is  unfavorable  to  the  slip  systems.  The  slips 
of  (10  0),  (01  0),  and  (001)  crystal  planes  become  more  difficult 
and  the  shear  modulus  (G)  thus  increases  to  89.57  GPa.  In  contrast, 
the  C\2  value  decreases  as  the  16(c)  sites  are  fully  occupied  by  Li+, 
and  this  behavior  not  only  leads  to  a  smaller  bulk  modules  but  also 
results  in  a  larger  (Cn-Ci2)  value.  Therefore,  both  the  plasticity 
and  ductility  of  the  Li2Ti204  become  poor  as  the  Li+  intercalation 
process  is  further  considered. 

The  bonding  of  systems  also  significantly  affects  the  thermo¬ 
dynamic  properties.  To  obtain  the  formation  enthalpy,  the  molar 
reaction  enthalpy  (A rHm)  is  derived  from  the  crystal  energies  first. 
Relying  on  the  thermodynamic  cycle  described  below 

Li(s)  +  2Ti02(s,  rutile)  ->  LiTi204(s)  (9) 

A rHm  can  be  calculated  by, 

A  rHm  =  AfHm(LiTi204)-AfHm(Li)-2AfHm(Ti02(rutile)) 

=  £(LiTi204)  -  £(Li)  -  2£(Ti02(rutile))  (10) 

where  E(Li),  £(Ti02)  (s,  rutile),  and  £(LiT204)  represent 
the  total  energies  of  lithium  (-190.18902954715  eV),  tita¬ 
nium  oxide  (-2483.2406419855  eV),  and  LiT204  crystals 
(-5158.56640998375  eV).  According  to  the  DFT  algorithm, 


the  molar  reaction  enthalpy  (A rHm)  is  calculated  to  be 
-182.723  kjmol-1.  Because  the  absolute  value  of  formation 
enthalpy  is  unable  to  determine,  some  experimental  data  is 
needed  to  derive  the  formation  enthalpy  for  LiTi204.  According 
to  the  thermodynamic  data  for  substances,  the  AfHm  value  for 
Ti02  (s,  rutile)  is  -944.0 ±0.8 kjmol-1  [25].  On  the  basis  of  the 
following  equation: 

AfHm(LiTi204)  =  ArHm  ±  AfHm(Li)  ±  2AfHm(Ti02,  s,  rutile) 

(11) 

The  formation  enthalpy  of  LiTi204  (A fHm)  is  finally  determined  to 
be  -2070.723  ±  1.6  kj  mol-1.  The  formation  enthalpy  of  LiT204  is 
obviously  lower  than  that  of  LiCo02  (-142.54 ±  1.69 kjmol-1 
[26]),  LiNi02  (-56.21  ±  1.53  kj  mol-1  [26])  and  LiMn204 
(-1380.9  ±2.2  kjmol-1  [27]).  Hence,  LiT204  has  a  higher 
thermodynamic  stability  than  that  of  LiCo02,  LiNi02  and  LiMn204. 

The  thermodynamic  cycle  of  Li2T204  described  below 

2Li(s)  ±  2Ti02(s,  rutile)  ->  Li2Ti204(s)  (12) 

A rHm  can  be  calculated  by, 

A  rHm  =  £(Li2Ti204)  —  2E(Li)  -  2£(Ti02(rutile))  (13) 

where  E(Li),  E(Ti02)  (s,  rutile),  and  E(LiT204)  represent  the 
total  energies  of  lithium,  titanium  oxide,  and  Li2T204  crystals 
(-5349.6273012625 eV).  The  molar  reaction  enthalpy  (ArHm)  is 
calculated  to  be  -267.080 kjmol-1.  On  the  basis  of  the  following 
equation: 

AfHm(LiTi204)  =  ArHm  ±  2Z\fHm(Li)  ±  2AfHm(Ti02,  s,  rutile) 

(14) 

The  formation  enthalpy  of  LiTi204  (AfHm)  is  finally  determined 
to  be  -2155.08  ±  1.6  kj  mol-1.  Hence,  it  can  be  concluded  that  the 
thermodynamic  stability  of  Li2Ti204  can  be  further  improved  after 


Fig.  2.  Electron  density  difference  diagram  for  LiTi204:  (a)  view  for  (001)  crystal  plane  and  (b)  three-dimensional  view. 
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the  Li+  intercalation  of  LiTi204.  The  high  thermodynamic  stabil¬ 
ity  and  good  ductility  of  LiTi204  are  expected  to  show  a  potential 
commercial  application  on  power  lithium-ion  battery. 

4.  Conclusions 

In  present  paper,  we  have  carried  out  systematic  investigation 
of  the  mechanical  and  thermodynamic  stability  of  spinel  LiTi204 
anode  material.  The  three  elastic  constants  (Cn,  Cj2  and  C44)  and 
mechanical  parameters,  such  as  shear  modulus  G,  bulk  modulus 
B,  and  Young’s  modulus  E  of  LiTi204  are  calculated.  Furthermore, 
the  effect  of  electronic  structures  on  thermodynamic  and  mechan¬ 
ical  properties  is  also  addressed.  The  calculated  results  suggest  that 
the  cubic  compound  LiTi204  is  mechanically  stable.  The  calculated 
G/B  ratio  of  0.584  indicates  a  good  ductility  of  the  LiTi204,  while 
(Cn-Ci2)  suggests  the  plasticity  of  the  LiTi204  is  relatively  poor. 
The  formation  enthalpy  of  LiTi204  shows  that  it  has  a  good  ther¬ 
modynamic  stability.  In  Li2Ti204  compound,  the  full  occupation  of 
16(c)  sites  by  Li+  leads  to  a  poor  plasticity  and  ductility,  while  the 
thermodynamic  stability  of  Li2Ti204  can  be  further  improved  after 
the  Li+  intercalation  of  LiTi204.  The  new  approach  presented  in  this 
work  is  promising  for  predicting  the  stability  of  electrode  materials. 
It  can  be  concluded  that  the  first-principles  approach  can  be  useful 
to  better  understand  the  stability  mechanisms  of  electrode  mate¬ 
rials  and  it  will  contribute  to  the  development  of  safe  and  stable 
electrode  materials. 
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